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We report a polarized neutron scattering study of the orbital-like magnetic order in strongly 
underdoped YBa2Cu3 06.45 and YBa2(Cuo.98Zno. 02)306.6- Their hole doping levels are located 
on both sides of the critical doping pmi of a metal- insulator transition inferred from transport 
measurements. Our study reveals a drop down of the orbital-like order slightly below pmi with 
a steep decrease of both the ordering temperature Tmag and the ordered moment. Above pmi, 
substitution of quantum impurities does not change Tmag, whereas it lowers significantly the bulk 
ordered moment. The modifications of the orbital- like magnetic order are interpreted in terms of 
a competition with electronic liquid crystal phases around pmi- This competition gives rise to a 
mixed magnetic state in YBa2Cu3 06.45 and a phase separation in YBa2(Cuo.98Zno. 02)306.6- 

PACS numbers: PACS numbers: 74.25-Ha 74.72.Bk, 25.40.Fq 
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There has been a long-standing debate among 
condensed-matter physicists about the origin of the 
pseudo-gap (PG) state in high-temperature supercon- 
ducting (SC) cupratesi. Two theoretical approaches have 
been opposed: in the former, the pseudogap state is a 
precursor of the SC state, in the later, its corresponds to 
another state of matter, competing with superconductiv- 
ity. Supporting the second scenario, polarized neutron 
diffraction has recently revealed the existence of a 3D 
long range magnetic phase, hidden inside the pseudo- 
gap state of underdoped YBa2Cu306-i-,;^'^. This phase 
is likely to be a generic feature of the pseudo-gap state, 
as it has been later observed in another cuprate fam- 
ily HgBa2Cu04+^. In both cuprates families, its or- 
dering temperature T^ag decreases linearly as a func- 
tion of hole doping (p) and vanishes near ppg^O.19 
(Fig. [T|) , the end point of the pseudo-gap phase accord- 
ing to thermodynamic measurementsS,. T^ag matches 
the pseudo-gap temperature T* deduced from resistiv- 
ity measurements'*'^ or appears slightly lower than T* in 
other resistivity-^ and NMR Knight shift studies* depend- 
ing of the criterion to determine T*. Thus, T* could be 
viewed as the onset temperature associated with the mag- 
netic ordering. In addition, a singular point occuring at 
Tmag was identified in susceptibility measurements car- 
ried out in underdoped YBa2Cu306-i-x and interpreted as 
a thermodynamic indication for the existence of a phase 
transition in the pseudogap stated. 

The magnetic phase can be described as a Q=0 an- 
tiferromagnetic (Q=0 AF) state: time reversal symme- 
try is broken, but lattice translation invariance is pre- 
served, since the same AF pattern develops within each 
unit cell. The occurence of such a magnetic state has 
been predicted in the circulating current (CC) theory 
of the pseudo-gap proposed by C.M.Varma^-. In this 
theory, staggered current loops give rise to orbital-like 
magnetic moments within Cu02 plaquettes. It is worth 



noticing that the staggered magnetization also couples 
to the uniform magnetization, yielding the weak sin- 
gularity in the uniform susceptibility at the ordering 
temperature-^, observed experimentally^. The overall 
symmetry of the magnetic neutron scattering pattern is 
actually consistent with the fourfold degenerated CC-6'// 
phase of Refj^^. 

Owing to the confinement of the current loops within 
Cu02 planes, orbital-like magnetic moments are ex- 
pected to be parallel to the c-axis. However, the ex- 
perimental observation of a substantial planar magnetic 
component^"— casts some doubt on the validity of the CC- 
phase as the origin of the Q=0 AF state. Nevertheless, 
several theoretical studies propose different ways to over- 
come this difficulty. For YBa2Cu306-i-x, for instance, a 
spin-orbit coupling effect can explain the appearance of 
a weak planar magnetic component^. In addition, an 
extremely weak ferromagnetism may accompany the CC 
order in magneto-optical measurements (Kerr effect 
For a CuOg octahedron (realized in the HgBa2Cu04+5 
system), variational Monte-Carlo numerical simulations 
show that circulating currents could be delocalised over 
all oxygen orbitals, even the apical onesi^. In the CC-6'// 
phase, there are 4 states corresponding to different config- 
urations of the current loops within the Cu02 plaquettes. 
In the initial CC-theory, this leads to 4 classical domains. 
Once excitations are taken into accounti^, it turns out 
that the ground state is made of the quantum super- 
position of the 4 states. Then, these quantum correc- 
tions can account for the observation of orbital-like mag- 
netism with a planar magnetic component, as observed in 
polarized neutron experiments. In addition, the theory 
predicts specific quasi non-dispersive collective magnetic 
excitations and one of them has been recently observed 
in polarized inelastic neutron scattering measurement in 
IIgBa2Cu04+a*^. Thus, to date, the CC-phase remains 
the only available theoretical approach able to account 
for the Q=0 AF order observed by polarized neutron 
diffraction. 
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FIG. 1: (color online) Schematic phase diagram of 
YBa2Cu306+x as a function of hole doping p; At low doping, 
the system is a Mott insulator and becomes an antiferromag- 
net (AF) below the Neel temperature Tjv (green area}^i^. 
Upon increasing hole doping, it enters the superconducting 
(SC) state below the critical temperature Tc (black bold line). 
p is deduced from ref,^^. The hatched superconducting area 
illustrates the reduction of the superconducting critical tem- 
perature Tc in YBa2(Cuo.98Zno.o2)306+x. The hole doping 
of both samples, that are studied here, are indicated in the 
phase diagram by a black point for Y645 (p=0.08, Tc=35 K) 
and a red point for Y66-Zn (p=0.12, Tc=30 K). 

While the Q=0 AF order appears 3D and is of long- 
range character in YBa2Cu306+5^^>^ and HgBa2Cu04-i-5*, 
its stability in La2-xSrxCu04 seems more questionable. 
In this system, for a hole doping p=0.085, the Q=0 AF 
order remains 2D and is short range^^. It further sets 
in below Tmag ^ 120 K, a temperature much smaller 
that the expected value of T* for such a low hole dop- 
ing. In this cuprate family, the expansion of the Q=0 
AF order is likely to be limited by the occurence of a 
competing electronic instabilityii. In La2-xSrxCu04, 
this competing phase could take the form of fluctuat- 
ing stripes which tend to become static close to a hole 
doping of p=l/8 or in presence of disorder (impurities, 
random potential, point defects) and external perturba- 
tions such as a uniform magnetic field— li^. Following 
numerical simulations22i2i indicating that a CC-phase 
could develop in two-leg ladders at least at short range, it 
has been suggested that currents loops could be confined 
within charge stripes in La2-xSrxCu04^. 

This observation rises questions about the stability of 
the orbital-like magnetic order in strongly underdoped 
YBa2Cu306+x- The order has been observed down to 
a hole doping of p=0.09. At slightly lower hole doping 
Pmi ^^0.085, a metal-insulator (MI) transition is inferred 
from quasiparticule heat transport^^ : this critical doping 
is assumed to be related to the freezing of stripe corre- 
lations at low temperature, as revealed by the in-plane 
resistivity anisotropy (the blue dashed area in Fig. [1]) . In 
addition, the transport data has been interpreted as a 
critical divergence of the cyclotron mass (to*), yielding 
a collapse of the Fermi temperature (oc 1/to*) around 



Pmt^ ■ However, a different conclusion has been recently 
drawn from other transport measurements^i where the 
metal-insulator crossover at pMi is claimed to be driven 
by a Lifshitz transition, namely it would correspond to 
a change in Fermi-surface topology at the critical con- 
centration Pa// where the electron Fermi pockets vanish. 
Without going into that specific debate, it is here im- 
portant to stress that all transport measurements point 
towards a critical doping, pmi, with different electronic 
properties on both sides. At the same doping level, it is 
worth to remind that the spin excitation spectrum also 
changes drastically^^"— i^. Likewise, the spontaneous 
appearance of an a-b anisotropy in spin correlations pro- 
vides strong evidences of electronic liquid crystal (ELC) 
phases near p}^J]2Ir^, Under the generic name of ELC 
phases, one usually group all states of correlated quan- 
tum electronic systems that break spontaneously either 
rotational invariance or translation invariance"^^. Using 
such a terminology, a fluctuating stripes state breaking 
the C4 rotational invariance would be classified as a ne- 
matic ELC phase, whereas a static stripes state that fur- 
ther break the translation invariance would correspond 
to a smectic ELC state. 

We here report a study of the stability of the 3D 
Q=0 AF order in YBa2Cu306+x for two distinct cases: 
(i) at low hole doping for p=0.08 [p < pmi), (ii) 
at larger hole doping p—0.12 (p > pmi) in the pres- 
ence of a disorder introduced through substitution of 
non magnetic Zn impurities (Fig [1]). Here we report 
a polarized neutron scattering study of the Q=0 AF 
state in two samples YBa2Cu3 6.45 (Y645, p=0.08) 
and YBa2(Cuo.98Zno. 02)306. 6 (Y66-Zn, p=0.12), where 
the fingerprints of an ELC phase has already been ob- 
tained in previous unpolarized inelastic neutron scatter- 
ing experiments^!^. 



II. EXPERIMENTAL DETAILS 

Sample and experimental set-up 

The samples are arrays of co-aligned tiny single crys- 
tals glued on a Silicon plate (Y645) or on an Al-grid 
(Y66-Zn). They are detwinned, with a detwinning ra- 
tion of 1/10 for Y645S2 and 1/4 for Y66-Zn29. Both 
samples consist of ~ 100 co-aligned and detwinned single 
crystals. For the Y645 sample array, the sample prepa- 
ration has been previously detailled in the supporting 
Online Material of ref,2i. was measured using SQUID 
magnetometry in a large fraction of the individual single 
crystals, yielding the average values of Tc—i^ K (Y645)— 
and Tc=30 K (Y66-Zn>22. Fig. [Ushows SQUID measure- 
ments or three representative single crystals. All crystals 
were characterized by magnetometry and were found to 
exhibit superconducting transition temperatures = 35 
K with transition widths of 2-3 K, testifying to their 
high quality. The bulk Tc of the mosaic sample was fur- 
ther crosschecked using the neutron spin depolarization 
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FIG. 2: (color online) SQUID susceptibility from 3 different 
representative single crystals constituting our sample array of 
YBa2Cu306.45 (p=0.08, i;=35 K). 



techniqueSfi on the (004) Bragg reflection. The lattice pa- 
rameters, a = 3.8388 A, b = 3.8747 A, and c = 11.761 A, 
correspond to the chemical composition YBa2Cu3 06.45, 
from which a hole concentration of 0.08 ± 0.005 per Cu02 
plane is deduced using Ref<^ and the Tc transition widths 
(Fig. [2|). Both samples exhibit thus a similar Tc. For 
Y66-Zn, the substitution of 2% Zn/Cu reduces Tc from 
64 K in pure Y66^ down to 30 K in Y66-Zn, while keep- 
ing the hole doping unchanged (p~ 0.12). For Y645, the 
small Tc value of 35 K is due to the reduction of the 
number of doped hole per Cu02 planes (p~ 0.08). 

Polarized neutron scattering measurements were per- 
formed on the triple-axis spectrometer 4F1 at the reactor 
Orphee in Saclay (France). The polarized neutron scat- 
tering set-up is similar to the one used in previous experi- 
ments on the same topic?"*^-: the incident neutron beam is 
polarized using a polarizing super-mirror and the polar- 
ization of the scattered beam is analysed using a Heusler 
analyzer. Standard XYZ-Helmholtz coils guide the neu- 
tron spin polarisation on the sample. The experimental 
set-up further includes, on the incoming neutron beam, a 
Mezei flipper for flipping the neutron spin direction, and 
a pyrolytic graphite filter for eliminating higher harmon- 
ics. For the polarized diffraction measurements, the inci- 
dent and final neutron wave vectors are set to 2.57 A~^. 
Samples were attached to the cold head of a 4K-closed 
cycle refrigerator and aligned in the [010]/ [001] scatter- 
ing plane, so that tranferred wave vectors Q of the form 
(0,K,L) are accessible. It is important to notice that, at 
variance with studies carried with twinned samples, the 
[100] and [010] directions are not equivalent anymore in 
the case of the present study on detwinned samples. 



Polarized neutron measurements and data analysis 

The magnetic signal we are looking for is induced by 
a Q=0 AF order, implying that the magnetic scattered 
intensity is superimposed onto the nuclear scattering on 
Bragg refiections. The magnetic intensity is typically 4 
orders of magnitude weaker than the strongest nuclear 
Bragg reflections. Therefore, polarized neutron tech- 
nique is essential to disentangle nuclear and magnetic 
contributions on Bragg reflections. Owing to the weak- 
ness of the magnetic intensity, the success of this mea- 
surement relies on the quality of polarization of the neu- 
tron beam, given by the flipping ratio (i?o)- To this end, 
the scattered intensity on a selected Bragg reflection is 
systematically measured in the spin-flip (SF) and non- 
spin-flip (NSF) channels, with 3 orthogonal neutron spin 
polarizations: Hx,Hy,Hz. For Hx and Hy, the neutron 
spin polarization is respectively parallel and perpendic- 
ular to Q in the scattering plane. For Hz, the neutron 
spin polarization is perpendicular to the scattering plane. 
In the rest of the manuscript, the indices SF,NSF and 
a=x,y,z indicate to which channel and to what kind of 
polarizations the measurements correspond. 

The magnetic intensity, Imag and the inverse flipping 
ratio R~^{T) are given by the following equation: 

K\T) = = i!!l^^ + i?-i(r) (1) 

^NSF,a ^NSF.a 

I stands for the Bragg intensity and Rq]^{T) is the inverse 
flipping ratio in the absence of any magnetic signal. 

The study of the magnetic signal is ususally performed 
on Bragg reflections (0,1, L) and (1,0, L) with integer L 
values. For twinned samples, (0,1, L) and (1,0, L) nuclear 
Bragg reflections are equivalent. At variance, using a 
detwinned sample, the nuclear scattering at (0,1, L) is 
significantly weaker than at (1,0, L), since the nuclear re- 
sponse associated with CuO chains is absent. In the SF 
channel, the polarization leakage from the NSF channel 
is therefore substentially reduced for (0,1,L), allowing a 
better detection of the magnetic scattering. 

According to theory, Rq ]^ should not depend on Q, 
neutron polarization and temperature. In practice, that 
is not the case due the imperfect experimental setup and 
a shift of the sample position within the guide field when 
changing temperature. Thus, for each Bragg reflection 
(0,1, L) and for a given polarization a , one calibrates 
the magnitude of Rq]^ at high temperature where the 
magnetic signal is absent. The temperature dependence 
of Rq ]^ is measured at the larger wave vector Q= (0,2,0) 
or (0,2,1) where the magnetic scattering is vanishingly 
small. 

In the SF channel, the scattered intensity is obviously 
dominated by this polarization leakage at all tempera- 
tures. Nevertheless, one may also need to pay attention 
for the background on top of which the Bragg scattering 
develops. Indeed, the temperature dependence of this 
background can blur the extraction of the temperature 
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dependence of the magnetic intensity. For each studied 
Bragg reflection (0,1, L), a systematic measurement of the 
background is therefore performed at (0,0.9,L) away from 
the Bragg scattering in both the SF and NSF channels 
and substracted from the scattered intensity measured at 
the Bragg reflection. This background substraction was 
not carried out in previous studies^''^: it actually corre- 
sponds to a refinement of the method used to extract 
the magnetic signal. It allows a better determination of 
the temperature dependence of the magnetic signal, es- 
pecially when this signal is weaker than in the previous 
studies: this is precisely the case in the present study. 

Once the magnetic intensity is determined, it can be 
calibrated in absolute units (mbarns) by comparison with 
the intensity of the (0,0,4) nuclear Bragg reflection (7 
barns): this intensity is indeed hardly dependent on the 
oxygen content and provides a good reference for calibra- 
tion. This calibration procedure is the same as the one 
used in previous studies^—. 

III. EXPERIMENTAL RESULTS 

Q=0 AF order at low hole doping 

To begin with, let us consider the hole doping effect on 
the stability of the 3D Q=0 AF order. Previous studies 
established that Tmag and the value of the order mag- 
netic moment at low temperature grew continuously in 
YBa2Cu306+x, when reducing the oxygen content from 
x=0.75 down to x=0.5^. In Y65 (Tc=54 K), the hole 
doping is p~0.09 (> pmi) and Tmag is close to room 
temperature^. The study of our Y645 sample allows us 
to track the evolution of the magnetic signal at lower 
oxygen content (x==0.45) and slightly lower hole doping 
(p=0.08). 

Figure |3]shows R^^{T) at the Bragg reflections (0,2,0), 
(0,1,0) and (0,1,1). At (0,2,0), R-^T) exhibits an al- 
most linear decrease on cooling down. This allows us 
to set the T-dependence of the reference inverse hipping 
ratio: R^l (x{l + AxT) with A = 5xlO-'^K-\ At vari- 
ance, at (0,1,0) and (0,1,1) , R~^{T) increases smoothly 
on top of R'al:{T) below T^ag ~ 170 K ± 30 K. Accord- 
ing to Eq. [TJ this indicates the gradual appearance of the 
magnetic intensity below Tmag- Figure [3] highlights two 
important features. For a small reduction of the hole 
doping, Tmag Steeply drops from 300 K (p > pmi) down 
to 170 ± 30K [p < Pmi)- Meanwhile, the magnetic inten- 
sity at (0,1,0) and (0,1,1) decreases from 9 ±1.8 mbarns 
and 2.9±0.3 mbarns in Y65^ down to 1.7±0.2 mbarns 
and 0.3±0.2 mbarns here in Y645. Furthermore, the T- 
dependence of the magnetic intensity suggests that Tmag 
should be viewed as an onset temperature rather than a 
net transition temperature. Our data therefore demon- 
strate that the 3D Q=0 AF order exhibits a sudden de- 
crease in strongly underdoped Y645, i.e slightly below 
Pmi, and vanishes at lower doping on approaching the 
insulating state. 
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FIG. 3: YBa2Cu306.45 - Temperature dependence of 

at different Bragg reflections. Solid lines are guides to the 

eye. The dotted lines indicate Rq]. oc (I + A x T) with A — 



It is generally believed that T* grows continuouly when 
decreasing the hole doping and should be much larger 
than room temperature upon approaching the Mott AF 
phase. In this case, our study shows that Tmag does 
not seem to match T* anymore at low doping. How- 
ever, the exact doping behavior of T* at such low doping 
might have been overlooked. It is indeed worth pointing 
out that the T* determination in strongly underdoped 
YBa2Cu306+x is not very accurate neither in NMR 
Knight shift§- nor with resistivity measurements^. In 
Rcf.^, samples with the nominal oxygen content x=0.45 
exhibit much higher Tc than our Y645, suggesting that 
their hole doping could be p r-^OA rather than p=0.08. 
Furthermore, measurements above room temperature are 
almost impossible owing to the onset of oxygen mobility 
in the samples. For resistivity measurements carried out 
on strongly underdoped samples^, the detection of an 
anomaly around our value of Tmag could be hidden by 
the upturn of the resistivity at low temperature which 
characterized these samples. Therefore, what we found 
here for Tmag might be the actual doping behavior for 
T*. 

For a doping level of p=0.085, the Q=0 AF order de- 
velops at a low temperature in the monolayer system 
La2_xSrxCu04 {Tmag=^'20 K)ii, while in bilayer system 
YBa2Cu306+x the order develops slowly Tmag ^170 K 
and is further enhanced below '^120 K (Fig. [3]). In the 
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former system, the order is 2D and at short range. In 
the later system, the magnetic neutron scattering inten- 
sity is still observed on Bragg reflections, suggesting that 
3D correlations survive: Attempts to detect a 2D mag- 
netic intensity at Q=(0,1,L) with non-interger L values 
do not reveal any signal. In La2-xSrxCu04, the body 
centered crystal structure, the intrinsic disorder due to 
Sr/La substitution and an antiparallel stacking of charge 
stripes along the c-axis could contribute to the reduction 
of the Q=0 AF correlations along the c-axis. 



Dilution of the Q=0 AF order through substttution of 
quantum impurities 

We study in this section the effect of non-magnetic 
Zn impurities on the Q=0 AF order for a doping level 
above pmi- While Tc is strongly reduced by non mag- 
netic impurity, neither hole doping nor T* are modi- 
fied. However spin correlations are known to be enhanced 
locally around Zn impurities, according to local probe 



(Oil) 



measurement 



,40,41 



The local spin correlations are ex- 
pected to compete with the Q=0 AF order. 

As in Y66'^, the T-dependence of the reference inverse 
flipping ratio, RqI.{T), is measured on the Bragg reflec- 
tion (0,2,1). Fig. Hla shows that RqI{T) is actually a 
constant for that particular experiment. Furthermore, 
R^^{T) is also temperature independent at Q=(0,0,2). 
This is a characteristic feature of the Q=0 AF order, 
which does not yield any magnetic scattering on Bragg 
reflections (0,0, L) with integer L values^'^^. At vari- 
ance, R^l{T) starts growing below rmQg=250±20K at 
Q=(0,1,1) (Fig. Sla) and (0,l,0)(Fig. Hb). At (0^1,0), 
one can check that the polarization sum rule R~^ = 
R~^ + R~^ (valid for a magnetic signal only in absence of 
any chirality) is fulflUed within error bars: This demon- 
strates of the magnetic nature of the signal that develops 

below Tmag- 

In Y66-Zn, we observe a magnetic signal reminiscent of 
the one reported in pure Y66^. It displays the same kind 
of structure factor, characterized by the appearance of a 
magnetic response on (0,1, L) Bragg reflections and the 
absence of a magnetic response at on Bragg reflections 
(0,0, L). The polarization analysis of the data is consistent 
with magnetic moments oriented at ~45° with respect to 
the c-axis, in agreement with similar analysis in the Zn 
free coumpound^i^. The temperature dependence can be 
described using the a power law (1 — T/Tmagf^ (with 
P = 0.185 ± 0.06) as reported for pure Y66 (Fig. H]). Fi- 
nally, the magnetic order sets in at almost the same tem- 
perature Tmag as in pure Y66^i^. This is in agreement 
with Knight NMR measurements, which demonstrated 
that the onset pseudo-gap temperature T* was not af- 
fected by substitution of non magnetic Zn impurities^. 

Nevertheless, dilute quantum impurities induce a net 
perturbation of the Q=0 AF order, which becomes ap- 
parent after a more quantitative comparison of data in 
Y66^ and Y66-Zn (Fig. The magnetic intensities at 
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FIG. 4: (color online) YBa2(Cuo.98Zno. 02)305. 6 - a) Temper- 
ature dependence of at different Bragg reflections. The 
dotted lines indicate Rq\{T) which is a constant for that set 
of measurements, b-d) Temperature dependence of R^^ at 
(0,1,0) for different polarizations: b) Hx (black), c) Hy(red), 
d) Hz (green). The solid lines indicate the appearance of the 
magnetic signal at Tmag ~ 250 K characterized by the tem- 
perature dependence oc (1 - T /Tmagf^ with /3=0.18^. For 
measurements on the Bragg reflections (0,2,1) and (0,0,2), the 
background is measured at Q=(0,2±0.1,l) and (0,0,2±0.2), 
respectively. 



low temperature on the (0,1,0) and (0,1,1) Bragg reflec- 
tions are 5.4±0.4 mbarns and 1.4±0.2 mbarns respec- 
tively in pure Y66. They drop down to 2±0.5 mbarns 
and 0.9±0.2 mbarns under Zn substitution. Since the 
neutron scattering measurement is a bulk measurement, 
the reduction of the magnetic intensity can be accounted 
for by a reduction of the sample fraction occupied by the 
magnetic order. In a dilution model, also called "swiss 
cheese model" the magnetic order within CUO2 planes 
is preserved far from Zn impurities: Figure [S] gives a car- 
toon picture of the inhomogenous magnetic distribution 
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FIG. 5: (color online) Array of 10x10 Cu02 plaquettes. The 
grid indicate the plaquettes where the Q=0 AF order should 
still be present. The black full dots stand for the Zn im- 
purities (3% per Cu in Cu02 planes) which are surrounded 
by an area where the Q=0 AF order is supposed to vanish 
(large orange disks). In addition, quasi- ID incomensurate 
spin fluctuations are restored at low energy around Zn im- 
purities: their characteristic correlation length ^ is about 2.9 
The yellow circles around impurities suggest a possible 
overlap between static Q=0 AF correlations and dynamical 
incomensurate spin correlations. 



induced by the development of protected clusters around 
local defects. Far from Zn impurities, the ordering tem- 
perature and the magnitude of the order moment remain 
those of the pure system: they are exclusively controled 
by the hole doping, which is not modified by Zn subti- 
tution. At variance, the magnetic order locally vanishes 
close to Zn impurities. Considering , on the one hand, 
an overall reduction of the magnetic intensity by a factor 
of 0.5±0.2 and , one the other hand, that the magnetic 
intensity in neutron diffraction measurements is propor- 
tional to the square of the ordered magnetic moment (M) , 
one can estimate that the magnetic order is present in 
70±10% of the Cu02 plaquettes. In a dilution model, 
the magnetic moment reads: 



M{z) = Af(0)(l -7z) 



(2) 



with z the number of impurities per Cu in Cu02 planes: 
in YBa2(Cui_yZny)306+x, z — i.e 3% in our case. 7 
stands for the efficiency coefficient which defines the spa- 
tial extension of the disorder introduced around a given 
impurity. We obtain a value 7=10±3, implying that a 
Zn induced disorder spreads up to the third nearest Cu 
neighbors around the impurity site (Fig. [5]), in perfect 
agreement with previous estimation from NM R^°i^^ and 
inelastic neutron scattering2S. 



IV. DISCUSSION 

Hole doping induced decay of the magnetic order versus 
dilution effect 

Fig.[5]reports the temperature dependence of the mag- 
netic intensity at (0,1,0) calibrated in absolute units in 
pure Y66, Y66-Zn and Y645: it summarizes our main 
observations. In Y66-Zn, T^ag is not affected by Zn sub- 
stitution and remains close to the value reported in pure 
Y66. The net reduction of the low temperature ordered 
magnetic moment can be ascribed to a reduction of ef- 
fective sample fraction occupied by the Q=:0 AF order 
owing to the local destruction of Q=0 AF correlation in 
the vicinity of quantum impurities. In this swiss cheese 
scenario, Tmag and the ordered magnetic moments are 
preserved far from Zn impurity, being essentially con- 
trolled by the hole doping level which is not modified by 
Zn substitution. In Y645, the low temperature order mo- 
ment is of the same magnitude as in Y66-Zn, but Tmag is 
strongly reduced. Imag(T) does not follow the power law 
(1 - T/T^ag?^ (with P = 0.185 ± 0.06) as observed at 
larger doping in Y66^ (Fig.H]). Imag increases more grad- 
ually on cooling down (Fig. [6]) . This suggests that there 
could be a distribution of ordering temperatures or that 
the Q=0 AF correlation length still increases below T„iag- 
In general, for magnetic systems, the low temperature 
magnetic moment M(oc ^Imag in neutron diffraction) — 
and the ordering temperature are expected be propor- 
tional. From Y65 to Y645, only ~20% of the magnetic 
intensity remains at low temperature, corresponding to 
a reduction of the low temperature ordered moment by a 
factor ~ 2.3. A scaling relationship M oc T^ag would im- 
plies a similar decrease of T^ag from 300 K (Y65) down to 
130 K (Y645). This value is consistent with the observed 
Traag — 170 ± 30 K in Y645. One can therefore conclude 
that the reduction of both T^ag and the low temperature 
ordered moment can be ascribed to a uniform weakening 
of the 3D Q=0 AF order below pm/=0.08. 



Competition with electronic liquid crystal phases 

In strongly underdoped La2-xSrxCu04 (p=0.085), a 
2D short range Q=0 AF order develops below T^ag ^ 
120 K^^. In addition, the evolution in the same tem- 
perature range of the incommensurate (IC) spin fluctua- 
tions around the planar AF wave vector Q^i^-=(0.5,0.5) 
is of particular interest. In twinned La2-xSrxCu04, 
these fluctuations can be observed at planar wave vec- 
tors Q/c= Qaf ±(<5,0)=Qaf ±(0,(5). They are usually 
interpreted as the flngerprint of fluctuation stripes. Sur- 
prizingly, the intensity of these IC fluctuations displays 
an upturn below T^ag and the IC parameter i5 starts in- 
creasing below Tmag like an order parameter. This study 
suggests that the vicinity of a spin and/or charge density 
wave state is the main cause of the limitation of Q=0 AF 
correlations. 
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FIG. 6; (color online) Temperature dependence of 
the magnetic intensity Imag at Q=(0,1,0) (or (1,0,0) for 
twinned samples): YBa2Cu306.45 (black - Y645, Tmag ~ 
170 K ), YBazCuaOs.e (red - Y66, T^^g ~ 235 K^), 
YBa2(Cuo.98Zno.o2)306.6 (black - Y645, T^ag ~ 250 K). 



Following this study, we can try to interpret the evolu- 
tion of the Q=0 AF order in Y645 in the light of the spin 
correlations around Qaf present in this material and 
measured on the same sampleS^:. In Y645, the low energy 
spin fluctuation spectrum is characterized by IC spin fluc- 
tuations that spontaneously develop a net a-b anisotropy 
below Telc ^ 150 K, at the planar wave vector Qic= 
Qaf ±(<5,0) ((5=0.045). At low temperature, these fluctu- 
ations freeze gradually, yielding a quasi-elastic magnetic 
signal. At the time scale of neutron scattering experi- 
ment, a quasi-lD short range spin density (SDW) order 
is detected below Tsdw ~ 30 K, but magnetic corre- 
lations become really static below a few Kelvin accord- 
ing to /iSR measurements^!. Telc is interpreted as the 
fingerprint of a nematic ELC phase, where the C4 rota- 
tional symmetry is spontaneously broken owing to strong 
electronic correlations^. At low temperature, the lat- 
tice translation is further broken in the SDW phase, that 
can be viewed as (glassy) smectic ELC phase. The com- 
plete study of spin correlations in strongly underdoped 
YBa2Cu306-(-x from x=0.3 to 0.45^, reveals a continuous 
decrease of Telc and Tsdw when increasing hole doping 
(Fig. [7]). The SDW phase should vanish when approach- 
ing the metal-insulator critical doping pM/ '^0.085, which 
seems to be the quantum critical point associated with 
the SDW state. Further evidence for quantum criticality 
is derived from the scaling properties of the dynamical 
spin correlations of Y645^. In line with this scenario, 
the scattering function at T = exhibits a spin gap for 
P > Pmi—1 but a dilute concentration of nonmagnetic Zn 
impurities locally restores the SDW in Y66-Zn2S. 

Figure [7] shows the evolution of the different tempera- 
tures Tmag, Telc ^ud TsDW as a function of hole doping 
for YBa2Cu3 06+x- When decreasing the hole doping in 
pure YBa2Cu306+x, Tmag decreases steeply across pmi- 




hole doping p 

FIG. 7; (color online) YBa2Cu306+x - Magnetic phase di- 
agram as a function of hole doping (p) (see Fig[T]). At low 
doping, when the system exhibits an insulating-like behav- 
ior at very low temperature, a nematic electronic liquid state 
shows up below Telc (ligth blue squares) and transforms 
into the spin density wave (SDW) state below Tsdw (dark 
blue triangles)-^''. In the metallic state, the Q=0 AF order 
appears below Tmag (black bullets)^'^. Dashed lines the lin- 
ear extrapolations of Tmag and Tsdw- The SDW may vanish 
around critical doping pm/ ~0.085. The Q=0 AF disappears 
at larger hole doping close to ppg ~0.19, which is also the end 
point of the pseudo-gap phase according to thermodynamic 
measurements^. The figure further shows that Tmag steeply 
decreases below Tmi- 



This indicates that the Q=0 order becomes destabilized 
by the ELC phase below pmi- In Y645 {p < pmi), the 
Q=0 AF order and the nematic ELC phase develop in the 
same temperature range: Tmag — 170 K and Telc '^150 
K. The appearance of static Q=0 AF correlation and 
the changes in IC spin fluctuations almost at the same 
temperature in both YBa2Cu06+x and La2-xSrxCu04 
around p ^0.08 suggest the existence of a generic in- 
terplay between the Q=0 AF order and the electronic 
instability generating the incommensurate spin correla- 
tions: it is conceivable that both phases coexist at low 
doping. In Y66 {p > pmi), the 3D Q=0 AF order is fully 
developped, but the ELC instability is still latent. As a 
result, a phase separation takes place in Y66-Zn when 
introducing local defects in Cu02 planes: the Q=0 AF 
order primarily vanishes around defects, while the SDW 
immediately nucleates around the Zn impurities (Fig.[S]). 
Far from Zn impurities, the Q=0 AF order remains es- 
sentially the same as in the Zn free system. Its ordering 
temperature and ordered magnetic moments are likely to 
be controled by the hole doping which remains unchanged 
under Zn substitution. 

At this stage, it is too early to give a definitive pic- 
ture of the mixed phase that could develop below pmi- 
Indeed, there is no clear consensus either about the ori- 
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gin of the Q=0 AF order, or about the origin of the 
quasi-lD IC-SDW phase. Up to now, the only model 
that could account for the Q=0 AF state is the circulat- 
ing current modcl^" with orbital-like magnetic moment 
expected primarily along the c-axis^^.. Concerning the 
SDW instability, one can notice that the evolutions of the 
incommensurate parameter 6 as a function of hole dop- 
ing are different in the monolayer system La2-xSrxCu04 
and the bilayer system YBa2Cu306+x^°- Such a differ- 
ence is not expected within a stripe model, but could 
be accounted for by the spiral-SDW model proposed by 
O.P. Sushko\«^. The orbital-like magnetic moments of a 
CC-phase and the Cu spins of planar spiral-SDW phase 
could mix owing to the spin-orbit coupling introduced in 
Ref.^^. As a result, Cu spins and orbital-like moments 
could be further tilted. In this scenario, a tilted spiral- 
SDW state would be expected. In addition, one can- 
not exclude a spatial modulation of the order parameter 
of the CC-phase. This scenario is still highly specula- 
tive and further theoretical and experimental works are 
needed to get a better understanding of such a mixed 
phase. 



strongly underdoped YBa2Cu3 06.45 and underdoped 
YBa2(Cuo.98Zno.o2)306-i-x- Our data indicate a dilution 
of the 3D Q=0 AF through Zn substitution, whereas, the 
Q=0 AF order undergoes a steep drop down when low- 
ering the hole doping below a critical hole doping level 
Pa//=0.085. Both phenomena suggest a competition with 
(nematic and smectic) electronic liquid crystal phases, 
previously observed in the same samples in inelastic and 
quasi-elastic neutron scattering measurements. While 
polarized neutron diffraction has highlighted the appear- 
ance of stagerred magnetic moments in Cu02 plaque- 
ttes, recent STM measurement in Bi2Sr2CaCu208+5— 
suggest that the electronic nematicity could also find its 
origin within the Cu02 plaquettes. Combining all the 
information should lead condensed matter physicists to 
reconsider the roles of copper and oxygen in the physics 
of cuprates and the way an effective single band model 
should be derived from the 3-band Hubbard model to 
account for the low energy charge and spin properties in 
these materials. 



CONCLUSION 



We report a polarized neutron diffraction study 
of the evolution of the 3D Q=0 AF order in 
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